Non-gaussianity from the trispectrum in general single field inflation 
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We compute the fourth order action in perturbation theory for scalar and second order tensor 
perturbations for a minimally coupled single field inflationary model, where the inflaton's lagrangian 
is a general function of the field's value and its kinetic energy. We obtain the fourth order action 
in two gauges, the comoving gauge and the uniform curvature gauge. Using the comoving gauge 
action we calculate the trispectrum at leading order in slow roll, finding agreement with a previously 
known result in the literature. We point out that in general to obtain the correct leading order 
trispectrum one cannot ignore second order tensor perturbations as previously done by others. The 
next-to-leading order corrections may become detectable depending on the shape and we provide 
the necessary formalism to calculate them. 



INTRODUCTION 



The theory of slow-roll inflation generically predicts that the observed cosmic microwave background radiation 
(CMBR) anisotropics are nearly scale invariant and very gaussian. Indeed, the latest observations of CMBR by 
WMAP3 [l| confirm these expectations. This constitutes one of the biggest achievements of modern cosmology. 

Despite its successes the theory of inflation still has many open questions. For example, we do not know the origin 
of the scalar field whose energy drives inflation, not to mention that we have never detected directly in the laboratory 
these kind of particles. The energy scale at which inflation happened is unknown by many orders of magnitude. 
There are many models of inflation that give similar predictions for the power spectrum of primordial perturbations, 
so which one (if any) is the correct one? 

For us to move a step forward in our understanding of the very early universe we have to work in two fronts. First 
the observational side. In the next few years, with improved experiments like the Planck satellite, we will measure 
the CMBR anisotropics to an incredible resolution. For example the observational bounds on the bispectrum (the 
three point correlation function of the primordial curvature perturbation Q will shrink from the present WMAP3 
value —50 < Jnl < 114 [l] to |/jvx| ^5 0, where the parameter Jnl parameterizes the size of the bispectrum. It's 
because this parameter is constrained to be small that we say that the CMBR anisotropics are very gaussian. The 
observational bounds on the trispectrum (four point function) will also tighten significantly from the rather weak 
present constraint of \tnl\ < 10 s [1,13] to the future constraint of |tjvl| ~ 560 ||, where tjvl denotes the size of the 
trispectrum. These previous observational bounds on the non-linearity parameters are for non-gaussianity of the local 
type. These bounds change depending on the shape of the wave vectors' configuration Q. This is one of the reasons 
why it is important to calculate the shape dependence of the non-gaussianity. 

In face of these expected observational advances, it is then imperative to push forward our theoretical knowledge of 
our theories and calculate more observational consequences of the different inflationary models to make a comparison 
with observations possible. One possible direction to be followed by us and many others is to calculate higher order 
statistics (like the trispectrum) of the primordial curvature perturbation. These higher order statistics contain much 
more information about the inflationary dynamics and if we observe them they will strongly constrain our models. 
Because these higher order statistics have a non-trivial momentum dependence (shape) they will help to discriminate 
between models that have a similar power spectrum (two point function). 

Calculations of the bispectrum for a single field inflationary model were done by Maldacena He showed that 
the primordial bispectrum is too small (of the order of the slow-roll parameters) to be observed even with Planck. 
Subsequent work generalized Maldacena's result to include more fields and more complicated kinetic terms [1, [^, [l(| • 
In Chen et al. have calculated the bispectrum for a quite general model of single field inflation. They showed 
that for some models even the next-to-leading order corrections in the slow-roll expansion may be observed. 

In this paper, we will focus our attention on the calculation of the trispectrum. In [ill ]. Seery et al. have calculated 
the trispectrum for slow-roll multi-field models (with standard kinetic terms) and they showed that at horizon crossing 
it is too small to be observed. But there are models of single field inflation, well motivated from more fundamental 
theories, that can produce a significant amount of non-gaussianity, such as Dirac-Born-Infeld (DBI) inflation [Hj[ll|]. 
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In [l3| the authors have computed the trispectrum for a model where the inflaton's lagrangian is a general function 
of the field's kinetic energy and the field's value, their result was obtained using a simple method [l4| , that only 
gives the correct leading order answer. In this paper we will provide the equations necessary to calculate the next-to- 
leading order corrections to the trispectrum. We argue that for some models these corrections might become equally 
observable in the future. In fact, we will calculate the fourth order action in the uniform curvature gauge that is exact 
in the slow-roll expansion and therefore in principle one could calculate all slow-roll corrections to the trispectrum of 
the field perturbations. 

We will also compute the exact fourth order action for the curvature perturbation £ in the comoving gauge. For a 
simpler inflation model (with the standard kinetic term) this was recently done in [l5[ . However [l5| did not consider 
second order tensor perturbations. We will argue that this is an oversimplification that leads to erroneous results. The 
reason simply being that at second order in perturbation theory, scalar degrees of freedom will source second order 
tensor perturbations and this will give a non-zero contribution for the fourth order action and hence the trispectrum. 

There are other reasons why we will perform the calculation in the comoving gauge. First of all, in doing so we 
work all the time with the gauge invariant variable C that is directly related with the observational quantities. The 
comoving gauge action can also be used other practical purposes. For example, it can be used to calculate loop 
effects that can possibly have important observational consequences. It can also be used to calculate the trispectrum 
of models where the potential has a "feature" (see [l|| [TtJ for an example of such calculation for the bispectrum). 
In the vicinity of the sudden potential "jump" the slow-roll approximation temporarily fails and one might get an 
enhancement of the trispectrum. There are well motivated models of brane inflation where the throat's warp factor 
suddenly jumps [HI]. 

This paper is organized as follows. In the next section, we introduce the model under consideration. In section ITTT1 
we shall study non-linear perturbations. First, we compute the fourth order action in the comoving gauge including 
both scalar and second order tensor degrees of freedom. After that we compute the fourth order action in the uniform 
curvature gauge. In section IIV1 we present the formalism needed to calculate the trispectrum. In section [Vj we 
calculate the leading order trispectrum using the comoving gauge action. We comment on previous works and on 
the observability of next-to-leading order corrections. Section IVT1 is devoted to conclusions. Finally in Appendix, we 
present the second order gauge transformation between the two previous gauges and a way to extract the transverse 
and traceless part of a tensor. 

II. THE MODEL 

In this work, we will consider a fairly general class of models described by the following action 

S=\f d^x^—g [M 2 Pl R + 2P(X, </>)] , (1) 

where <f> is the inflaton field, Mpi is the Planck mass that we will set to unity hereafter, R is the Ricci scalar and 
X = —^g^d^dvcf) is the inflaton's kinetic energy. t^ t „ is the metric tensor. We label the inflaton's lagrangian by P 
and we assume that it is a well behaved function of two variables, the inflaton field and its kinetic energy. 

This general field lagrangian includes as particular cases the common slow-roll inflation model, DBI-inflation 
P3] and K- inflation [2cj|. 

We are interested in flat, homogeneous and isotropic Friedmann-Robertson- Walker universes described by the line 
element 

ds 2 = -dt 2 + a 2 (t)S i3 dx l dx 3 , (2) 
where a(t) is the scale factor. The Friedmann equation and the continuity equation read 

3H 2 = E, (3) 

E = -3H (E + P), (4) 
where the Hubble rate is H = a/a, E is the energy of the inflaton and it is given by 

E = 2XP X - P, (5) 
where P x denotes the derivative of P with respect to X . 
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It was shown in [2l[ that for this model the speed of propagation of scalar perturbations ( "speed of sound" ) is c s 
given by 

2 = Px Px , fil 

We define the slow variation parameters, analogues of the slow-roll parameters, as: 

H_ = X£ 2 x m 
( H 2 H 2 ' 1 ' 

(9 » 

We should note that these slow variation parameters are more general than the usual slow-roll parameters and that 
the smallness of these parameters does not imply that the field in rolling slowly. We assume that the rate of change 
of the speed of sound is small (as described by s) but c s is otherwise free to change between zero and one. 

It is convenient to introduce the following parameters that describe the non-linear dependence of the lagrangian on 
the kinetic energy: 

Z = XP X + 2X 2 P, XX = (10) 
A = X 2 P XX + ^X 3 P XXX , (11) 

n-x 3 P™ + ^ 4 Px™. (12) 

5 

These parameters are related to the size of the bispectrum and trispectrum. The power spectrum of the primordial 
quantum fluctuation was first derived in [21[ and reads 

k ~ 36n 2 E + P ~ 8^ Cs£ ' (13) 
where it should be evaluated at the time of horizon crossing c ssr k = a*H*. The spectral index is 

dlnPf , , 

n.-l = -— L = -2e - v ~ 14 
rfmK 

WMAP observations of the perturbations in the CMBR tell us that the previous power spectrum is almost scale 
invariant therefore implying that the three slow variation parameters should be small at horizon crossing, roughly of 
order 10~ 2 . 

III. NON-LINEAR PERTURBATIONS 

In this section, we will consider perturbations of the background ^ beyond linear order. There is a vast literature on 
second order perturbations that are important when one is interested in calculating three point correlation functions, 
see for example @, H QU HH. In the present paper, we are interested in non-gaussianities that come from the 
trispectrum and so we need to use third order perturbation theory. For that we need to compute the fourth order 
in the perturbation action. In this section we will obtain the fourth order action in two different gauges. As a check 
on our calculations we will compute the leading order (in slow roll) trispectrum in both gauges. We will follow the 
pioneering approach developed by Maldacena 0] and used in several subsequent papers [a SEL G3- 
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For reasons that will become clear later it is convenient to use the ADM metric formalism (23[. The ADM line 
element reads 

ds 2 = -N 2 dt 2 + hij (dx l + N l dt) (dx j + N j dt) , (15) 

where N is the lapse function, N % is the shift vector and hij is the 3D metric. 
The action JT]) becomes 

S = 2 / dtd3x ^ N ( {3)R + 2P ) + \ J dt^xVhN- 1 {E l3 E ij - E 2 ) . (16) 
The tensor E^ is defined as 

= i (hij - ViNj - VjN^ , (17) 

and it is related to the extrinsic curvature by Kij = N~ 1 Eij. Vi is the covariant derivative with respect to hij and 
all contra-variant indices in this section are raised with hij unless stated otherwise. 
The hamiltonian and momentum constraints are respectively 

( '% + 2P -2n 2 N- 2 P^ x - N- 2 (E^E 13 - E 2 ) = 0, 

Vj (N~ 1 Epj — Vi (N~ 1 E) = nN^Vi^Px, (18) 

where 7r is defined as 

7T = <f>- N J Vj(f>. (19) 
We decompose the shift vector TV* into scalar and intrinsic vector parts as 

Ni = Ni + d t il), (20) 

where diN 1 = 0, here indices are raised with <5y. 

Before we consider perturbations around our background let us count the number of degrees of freedom (dof) that 
we have. There are five scalar functions, the field <j>, N, t/j, deth and hij ~ didjH, where H is a scalar function 
and deth denotes the determinant of the 3D metric. Also, there are two vector modes N l and hij ~ diXj, where 
X 3 is an arbitrary vector. Both iV* and x J satisfy a divergenceless condition and so carry four dof. Furthermore, 
we also have a transverse and traceless tensor mode 7^ that contains two additional dof. Because our theory is 
invariant under change of coordinates we can eliminate some of these dof. For instance, a spatial reparametrization 
like x l — x 1 + d l e(x,t) + e^(i,t), where e and eL are arbitrary and 9ieL = 0, can be chosen so that it removes one 
scalar dof and one vector mode. A time reparametrization would eliminate another scalar dof. Constraints in the 
action will eliminate further two scalar dof and a vector mode. In the end we are left with one scalar, zero vector and 
one tensor modes that correspond to three physical propagating dof. 

In the next subsection we shall use two different gauges that correctly parameterize these dof. Because physical 
observables are gauge invariant we know that both gauges have to give the same result for the trispectrum for instance. 
It seems then unnecessary to perform the calculation twice in different gauges. In practice, we will see that both gauges 
have advantages and disadvantages and one is more suitable for some applications than the other. Furthermore, it 
provides a good consistency check on the calculation. 



A. Non-linear perturbations in the comoving gauge 

In this subsection, we will compute the fourth order action for the general model ([T]) in the comoving gauge. In 
this gauge the scalar degree of freedom is the so-called curvature perturbation £ that is also gauge invariant. There 
are a few works on this subject using this gauge, see e.g. [l5j |. where the authors have calculated the fourth order 
action for a standard kinetic term inflation but they neglected second order tensor perturbations. We will show that 
this is an oversimplification that may lead to an erroneous result for the four point correlation function. 

In the comoving gauge, the inflaton fluctuations vanish and the 3D metric is perturbed as 



5(j> = 0, 

hij = a 2 e 2C -hij, h tj = 8 lj + 7^ + Ijjikrf + 



(21) 
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where deth = 1, 7^ is a tensor perturbation that we assume to be a second order quantity, i. e. 7^ = 0(£ 2 )- It obeys 
the traceless and transverse conditions 7? = 9 l 7y = (indices are raised with 6ij). C 1S the gauge invariant scalar 
perturbation. In (j2Tj) . we have ignored the first order tensor perturbations ^lijQw This is because any correlation 
function involving this tensor mode will be smaller than a correlation function involving only scalars, see results of 
Q- In the literature the second order tensor perturbations are often neglected, however they should be taken into 
account. The reason for this is because at second order the scalars will source the tensor perturbations equation. 
Later in this section, we will elaborate further on this point. Higher order tensor perturbations, like ^7ij, do not 
contribute to the fourth order action. 

We expand N and N l in power of the perturbation £ 



N = 1 + ai + a 2 + 



(22) 

(23) 
(24) 



where a n , Ni and ijj n are of order £ n . 

Some useful expressions for the quantities appearing in (|18[) . valid to all orders in perturbations but for 7^ = 0: 



= -2a- 2 e- 2C (dtCdX + 2d i d i () , 



(25) 



EijE ij -E 2 = ~6[H + ( 



b + d k (d k ij + d k (N k 



\m 3 (l 



d % N 3 + 8 3 N l + 28 t 8 3 \pd l N + d l 8^d l 8 3 ^ - d 2 ipd 2 ip 



-2diN 3 (d j (N l + d l (N j ) - 4N k d i Cd i d k ip - 20^ (d 3 (d 1 ^ + d l (d 3 4>) 
U),i),ci) l GYc + 2d 3 ( (djCNiN* + 2d j C,N i d i i) + djCW&i/f) 



(26) 



- 2 (h + () d,N + ■ 



-TV" 



-AT 



— (~d j N (diNj + djNi) + 2djN (d^Ni + d l QN 3 ) + 2diNd 2 ij) 

- 2d,C + a- 2 e- 2 < ( i (tfCdjNi - <)X<>,\ *) + ,U7>\'.V, - <U7>\\\, + did^Nj + d 2 ^ - ^d 2 N. t 



(27) 



In the previous equations, indices in the right-hand side are raised with while indices in the left-hand side are 
raised with hij . In the rest of this section indices will be raised with 5ij . 

Now, the strategy is to solve the constraint equations for the lapse function and shift vector in terms of £ and then 
plug in the solutions in the expanded action up to fourth order. 

At first order in £, a particular solution for equations (fl8|) is 



£ ;\7 (1) n ; 
a\ — — , Ni = 0, ipi 



C o2 2 £ }■ 



(28) 



At second order, the constraint equation for the lapse function gives 

—d 2 ip 2 = -2a- 2 d i ((dX + 2Hd l ip 1 ) - 4a x (a~ 2 ( - 2^() - 2a 2 (£ + 6A) 
-a- 4 {d i dkMd k il)i - d 2 ilj 1 d 2 ip 1 ) + 4a 2 (S - 3H 2 ) , 
and the equation for the shift vector gives 

2Hd ia2 - ^a- 2 8 2 N? ] = -a- 2 (d k a x d k d^i - d^d 2 ^ + d\d^i + 8^(8^) 



(29) 



(30) 



6 



Due to the fact that N l is divergenceless and that any vector can be separated into a incompressible and irrotational 

~(2) 

part one can separate the contributions from a 2 and N% in the previous equation. The irrotational part of Eq. (1301) 
gives 



2Ha 2 = 9" 2 9 l F l , 



(31) 



and the incompressible part gives 



=-d- 2 F i + d-%d k F k , 



(32) 



where Fi is define as the right-hand side of equation (|30|) . The operator d~ 2 is defined by d~ 2 (d 2 <f) = (p and in 
Fourier space it just bring in a factor of — 1/fc 2 . 

It was shown in [l(| that to compute the effective action of order ( n , within the ADM formalism, one only needs to 
use the solution for the Lagrange multipliers N and N l up to order ( n ~ 2 . Therefore in order to calculate the fourth 
order effective action the knowledge of the Lagrange multipliers up to second order is required. It is given in equations 
(25). (51) and J32H- 

The second order action is 



= / dtd 3 x 



a^e-ae{dCY 



The third order action is 0, H 

S3 = / dtd 3 x 



3a 3 e 



cc 2 



Za \ ti I a 



(33) 



(34) 



The scalar fourth order action is 



- / dtd 3 xa 3 
2 



-a~ 2 eC 2 (d() 2 + aj ( 2E + 9A + yh ) - 6Ca? (E 2A) 'K'n'pZ 



-2a 2 (E - 3H 2 ) + a- 4 ^ + C«i + a 2 J (d k d k N^ - (7'hf ) k .\f ) 

-2a- 4 (C + ax) d k Ni {1) (%ivj 2) - S kj 8 n N^ - 28^^ 
+ a - i (-45 ife iVj 1) S J 'CiV fe(2) + 2N k :i) d j (N k(1) d j () 



(35) 



Here, no slow-roll approximation has been made. The previous action has to be supplemented with the action 
containing terms with one and two tensors 



5 7 2 = ~ J dtd 3 x [a%jj ij - ad kllJ d k ^] , 



(36) 



S iC 2 = 



dtd 3 x 



(37) 



1. Canonical variable for quantization 



In order to calculate the quantum four point correlation function we follow the standard procedure in quantum 
field theory. However there is an important subtlety here. The gauge invariant quantity £ is not the correct variable 
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to be quantized, because it is not a canonical field . The canonical field to be quantized is the field perturbation Scj), 
or a convenient parameterization Cn defined by 

Cn = -£-8<t>, (38) 

where <fio is the background value of the field. We will see that C is related to Cn by a non-linear transformation so 
for the power spectrum calculation both procedures of quantizing C or Cn give the same answer because the difference 
between these two variables is a second order quantity. However, for the calculation of higher order correlation 
functions (like the bispectrum or trispectrum) Cn is the correct variable to be quantized as it is linear in 6<f>. For 
the trispectrum, quantizing C gives different results from quantizing Cn- We will find the relation between C and 
Cn through the gauge transformation equations from the uniform curvature gauge (discussed in detail in the next 
subsection) to the comoving gauge. In the uniform curvature gauge the ansatz is 



4>{x,t) = 4>a +S(j>(x,t), 



1 



hij = a 2 h i3 , hij = S i: j + % + - Jikjj H (39) 

where det/i = 1 and 7y is a tensor perturbation that we assume to be a second order quantity, i.e., 7y = 0(5(f> 2 ). 
The gauge transformations are 

C(x) = Cn(x) + J F 2 (C„(x))+^ 3 (C„(x)), (40) 
lij = 7*7 (*) + fMj, (41) 

where F 2 (^ n ),^ij = O(Cn), ^KCn) = C(Cn) are the terms coming from the second and third order gauge transforma- 
tions respectively, they can be found explicitly in Appendix [X] or in 0, HH • 

We now need to find the fourth order action for the variable Cn- Schematically, we write the different contributions 

as 

S Cn = S 4 (Cn) + MCn) + V 2 (Cn) + + S 2 (F 2 (C„)). (42) 

The first three terms come from Eqs. (|35H37p when we substitute C with Cn and 7,j with Eq. (|4"Tj) . Due to the 
non-linear relation between C and Cn, the third order action for C, Eq. (|34| . will after the change of variables give 
a contribution to the fourth order action like S3 (F2 (£„)). Similarly, the second order action, Eq. (j3U|) , will also 
contribute with S2(-p2(Cn))- I n principle, one would also need to compute the third order gauge transformation as 
the second order action gives origin to fourth order terms like Cn-FHCn) (where ^(Cn) is the third order piece of the 
gauge transformation). Fortunately, terms involving F3 can be shown to be proportional to the first order equations 
of motion for Cn, therefore when computing the trispectrum these terms will vanish and we don't need to calculate 
the third order gauge transformation explicitly at this point. It can be easily seen that equations (|34|) . ([35]) . (I36[) . 
(|3T|) or their counterparts in terms of Cn, Eq. (|4"2")) . have terms that are not slow roll suppressed. However, because 
in pure de Sitter space C is a gauge mode we expect the action (|4"2")l to be slow-roll suppressed (of order e). One can 
perform many integrations by parts to show that the unsuppressed terms of (|42p can be reduced to total derivative 
terms and slow-roll suppressed terms given by 



Mparts 



J dt d 3 x[ ~^[\ (djtn&Q 2 + l ~£d 2 (cKnC^Cn) + 9, Cr^" W (^Cn^Cn 



5e 



+ 23 4 (c^Cnd'Cn) d- 2 d t (dfcCnC^Cn) + \ (d^djdk (^Cn^Cn) ) ' 
dlCnd^nd"^ (dkCACn) - -/),.:,"\ h i) J (dfcCn^Cn) 

~d k dl (dkCnd.Cn) d- 2 d m l (d m CndlC n ) ] }■ (43) 



8a 3 H 4 



For us to be able to obtain the previous result it is crucial to include the contributions from the tensor actions 
and (|37[) . otherwise the trispectrum calculated using the O(e ) of (|33|) does not vanish, giving the wrong leading order 
result. Neglecting tensor perturbations (sourced by the scalars) for the calculation of the trispectrum is not consistent 
and leads to wrong results. This is one of the results of our work. The contribution from 7^ that comes through 7^- 
in (|36|) and (I37|) will result in terms that are already slow-roll suppressed and no further integrations by parts are 
need on these terms (see Eq. (15^1) of next subsection) . The final action for Cn is then 

S4£ n = S^ + Sparts, (44) 
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where denotes the terms of (|42j) that are suppressed by at least one slow-roll parameter. This final action is 
slow-roll suppressed as expected and no slow-roll approximation was made so it is also exact. 



B. Non-linear perturbations in the uniform curvature gauge 

In order to calculate the intrinsic four point correlation function of the field perturbation we need to compute the 
action of fourth order in the perturbations. In this subsection we will obtain the fourth order action in the uniform 
curvature gauge. In this gauge, the scalar degree of freedom is the inflaton field perturbation (50(x u ). There are 
several works in the literature where the authors also calculate the trispectrum. In [ll|, Seery et al. calculate the 
trispectrum of a multi-field inflation model, however the result is only valid for fields with standard kinetic energies, 
i.e. P(Xi, . . . , X n , 0i, ... , 0„) = X\ + ■ ■ ■ + X n — V, where X n is the kinetic energy of 0„ and V is the potential. 
In this paper we will generalize their result for an arbitrary function P(X, 0) but for single field only. Recently, 
Huang and Shiu have obtained the fourth order action for the model under consideration {!]). However the result was 
obtained by only perturbing the field lagrangian. This procedure gives the right result as long as we are interested in 
the leading order contribution in the small speed of sound limit and in the slow-roll approximation, which was their 
case of interest. In the present section we will compute the fourth order action that is valid to all orders in slow roll 
and in the sound of speed expansion. To do that we have to perturb the full action (JTJ) up to fourth order in the field 
perturbations. The procedure to obtain the fourth order action in this gauge is similar to the one used in subsection 

MM 

In this gauge, the inflaton perturbation does not vanish and the 3D metric takes the form 



(x,t) = 0o + S(j)(x,t), 

1. 
2 



h i3 = a 2 hij, h i:j = 8ij + 7-y 4- -7*7* H (45) 



where deth — 1 and 7^ is a tensor perturbation that we assume to be a second order quantity, i.e., 7^ = 0(<50 2 ). 
It obeys the traceless and transverse conditions 7* = d l %j — (indices are raised with 5ij). In the literature the 
second order tensor perturbations are often neglected, however based on our results it should be taken into account. 
We can always use the gauge freedom at second order to eliminate the trace and the vector perturbations of h%j ■ The 
presence of 7^ makes the three dimensional hypersurfaces non-flat so using the name uniform curvature gauge might 
be misleading. We will continue to use that name because in the literature that is the name given to the gauge where 
5(j)(x,t) = 0. 

We expand N and N l in powers of the perturbation <50(a;, t) 

N = 1 + ai + a 2 4 , (46) 

N, t = N t (1) + + ■ ■ ■ , (47) 

tp = + V>2 + • • • , (48) 

~ (™) 

where a n , Ni and ip n are of order (50™ and 0o(£) is the background value of the field. At first order in <50, a 



particular solution for equations (|T8|) is [7J, : 

n , = ^4> S4>Px, Ni (1) = 0, aVi = {-^H ) . (40) 



At second order, the constraint equation for the lapse function gives 
4 FT 1 

— d 2 ^ 2 = — (d 2 i>id 2 ^i - didj^d^xjji) + %H 2 (3a 2 - 2a 2 ) 



+ ^T ld2 ^ ~ (^o 2 al - 20 o 2 a 2 4- - 40 O (50a 1 ) P x - \ (d^d'S^ - 20 o 9^05Vi) P,x 
+<50 2 P 4- 20o(50 U ai - 6<p) P x<j) 

(-lOa 2 0o(50ai - 20 o d i (50<9 4 i/'i 4- 6a 2 o 2 a^ - 20 o 2 a 2 a 2 4- 4a 2 <50 2 - di508(^j P X x 
-0o (0o a\ + (50 - 20 O (50aiJ Pxxx 4- 20o <50 ( o ai - (50J Pxx^ - 0o (50 2 -P,x0</,, (50) 



and the equation for the shift vector gives 

1 -2c2aV( 2 ) 



P 



2Hd t a 2 - -a~ 2 d 2 M 2) = 4 ai Hd iai + a' 2 (d l a 1 d 2 ij 1 - A aiQjd fe Vi) + (S(j> - <j> ai) 

+di54>4>o (s'cf) - foai) Pxx + 8iS(j)4>oS(j)P^x4,, (51) 

where P^ means derivative of P with respect to <f>. 

Due to the fact that N l is divergenceless and that any vector can be separated into a incompressible and irrotational 

~( 2 ) . . . . 

part one can separate the contributions from a 2 and iVj in the previous equation. The irrotational part of Eq. [|5ip 



gives 

2Ha 2 = d~ 2 d l F u (52) 

and the incompressible part gives 

^2 Jvf ) = -0- 2 F t + D- A d t d k F k , (53) 



S A = / rf^ 



where F{ is defined as the right-hand side of equation ([51]) . 

The scalar fourth order action, where no slow-roll approximation has been made, is 

Si = S A + S B , (54) 

where 

- ^ [a 2 ai 5<p 2 + 2di50^ (S<t> - ax^o) + 2^d l 8cp (n? ] + ftVa) + a 2 ^ 2 a 3 + a^dS^P^ 

-— [4$a 4 a 2 (a 2 - 2a?) + 3</>oa 4 a 4 - 18$a 4 ar?<ty 2 + 80oa 2 <9 l c50 (ivf ) + 9^ 2 ) (V</> - <foai) 
-tyodiStfrd'ipi {[d8<j)) 2 ~ 3^a 2 a 2 + 8</> a 2 ai<^ - 3a 2 <50 2 ) + I2<p a 4 ai 5<j? + 4</> 3 1 a 4 a 3 50 
-4 2 a i (50aVic'fc(S0a fe Vi - 2</> 2 a 2 a 2 (<9£<?(>) 2 - 4a 1 ^<50a 2 (c«0) 2 - a 4 - a" 2 ^) 2 ) 2 ] P XX 

+ ^[6^ Q d t 5(t)d 1 ijj 1 (^foonStj) - Set) 2 - 0o«i) + 3(<9<S</>) 2 (2^1^ - ^a 2 ) - 160 o a 2 
+244> 2 a 2 a 2 5(j) 2 - 12$|a 2 a?50 + + 3a 2 ^ 2 (V/ - a - 2 (<9<S</>) 2 ) } P X xx 

+^^[a 2 5' ( f > 3 + ((9^) 2 + 24> o diS0ij^ (fa ai - S'ty + 30 2 a 2 a 2 <50 - 40 o a 2 a 1( ^ 2 ]Pxx0 

[2{d8<t>) 2 - 2a 2 8(t> 2 + 4a 2 ai0 o ^ + 2a 2 2 a 2 + ^o^^JPx^ 



+ ^o« 3 [60o« 2 <^ 2 ~ 40 o ai(5'</» 3 - 4^ga?«0 + + S<f> 4 ]P X xxx 
-jUf^a 3 (-54? + 30 o ai^ 2 + $a<? - 3$a?<ty) P,xxx</> 

+ ^(/) 2 ^ 2 a 3 ^-2ai^Q(5^ + 5(f) 2 + </> 2 a 2 ) P,xx<f>4, - ^4>aS(f> 3 a 3 (ai<fro - sty P,X4><s><s> , (55) 



Sb = / dtd 



/[ 1 1 1 

did 3 a: [afa^c^ + -a 3 a\54> 2 P^ + -a 3 ai<S</> 3 P 0<w , + — a 3 ^ 4 P 4 

^ [ - {di80i) X f - 2a l 5<j)a 2 d l 5<i)d l Tpi + 2a 2 cW (ivf ) + fl^) - «1<A)) 



2a 

+a 2 a 2 (d<50) 2 + 2 a 4 a 2 (a 2 - 2a 2 ) ] P x 



[29 l 7vi 2) %iV ?) (2) - A ai did k ^i [d i d k %l) 2 + e^N^\ + 12a 4 P 2 a 2 (a a - 2a 2 ) + Aa^^d 2 ^] . (56) 
The previous actions should be supplemented with the pure tensor terms and the tensor-scalar coupling terms: 

S = ~ J dtd 3 x [a%^ - ad k %d k ^\ , (57) 
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S = / dtd 



-fl ... 



(58) 



This constitutes the main result of this subsection. It is a good check for our calculation to see that the previous 
action (154|) reduces in some particular cases to previously know results present in the literature. 

For example, if we restrict our model to the standard inflation case, i.e., P(X, (j)) = X — V(4>), where V(tf>) is the 
inflaton potential, then all the terms in the scalar action (|55p vanish and the only contribution to the fourth order 
action comes from (|56|) . These terms exactly reproduce the result of Seery et al. [ll(, their equation (36), restricted 
to single field. However, in the total fourth order action there are also the tensor contributions (|57)l and (f58"|) . In 
general, to proceed one has to calculate the equation of motion for the second order tensor perturbations 7^ from 
eqs. (j5Il), dSll) to get 

a' / . \TT 

% + 2-% - a 2 7 y = {2P tX d j 6<l>a i 5(f> + iPxcPodjS^di^j , (59) 

where TT means the transverse and traceless parts of the expression inside the parenthesis (sec Appendix [B] for 
details of how to extract the TT parts of a tensor) and then solve this equation to obtain 7^ as a function of Sip. 
One can immediately see that at second order the scalars will source the tensor perturbation equation as it was 
previously shown by others 0, Hf| . At this order in perturbation theory, equation (|B"9"]) should also have a source 
term quadratic in the first order tensor perturbations, ^lij GW - We neglect these terms because we expect that any 
correlation function where ^lijQ W enters, which is sourced by the first order tensor modes squared, must be smaller 
than a correlation function with only scalars, see Ref. [?J for an example. In Fourier space, the source term of (f59"|) is 
suppressed by k 2 , where k is the wave number. Once we have the solution of 7^ in terms of 6<f) we can plug back the 
result in ([ST]) and {55]) to get the total fourth order scalar action. 



IV. THE GENERAL FORMALISM TO CALCULATE THE TRISPECTRUM 



A. The trispectrum of £ n 



Now we shall give the basic equations needed to calculate the trispectrum 0, [26| . First we need to solve the second 
order equation of motion for (obtained from (|33p h Defining new variables 



v k = zu k , z 



aV2e 



(60) 



where the Fourier mode function Uk is given by 
the equation of motion for £ ra is 



(61) 



// 1 2 7 2 * 
v k +c s k v k v k 

z 



0. 



(62) 



where prime denotes derivative with respect to conformal time r. This is also known as the Mukhanov equation. The 
previous equation can be solved, at leading order in slow roll and if the rate of change of the sound speed is small 
[101 ] . to give 



u k = u(r,k) 



iH 



y/Aec s k 3 



(1 + ikc s r) e 



—ikc s r 



(63) 



We do not need to impose any constraints in the sound speed and it can be arbitrary Only its rate of change is 
assumed to be small. The next-to-leading order corrections to the previous solutions are also known and can be found 
in [l(J. In the general case, we would have to solve Eq. (|6"2"]) without assuming slow roll. This can be done numerically. 

In order to calculate the £ n correlators we follow the standard procedure in quantum field theory. The curvature 
perturbation is promoted to an operator that can be expanded in terms of creation and annihilation operator as 



C„(t, k) = u(t, k)o(k) + u*(t, -k)at(_k). 



(64) 
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The standard commutation relation applies 

[a(kO, a f (k 2 )] = (2^) 3 <5( 3 )(k 1 - k 2 ). (65) 
The vacuum expectation value of the four point operator in the interaction picture (at first order) is 0, HfiJ 

(n|C„(t,k 1 )C„(t,k 2 )C„(t,k 3 )Cn(t,k4)|^) =-i f dt(0\ [Cn(t,k 1 )C„(t,k 2 )C„(t,k 3 )Cn(t,k4), J ff/(?)] |0), (66) 

Jto 

where io is some early time during inflation when the inflaton vacuum fluctuation is deep inside the horizon, t is 
some time after horizon exit. |fi) is the interacting vacuum which is different from the free theory vacuum |0). If one 
uses conformal time, it's a good approximation to perform the integration from — oo to because r rs —(aH) . Hi 
denotes the interaction hamiltonian and it is given by Hi = 7r£„ — L, where n is defined as n = J^r- and L is the 
lagrangian. In this work, we will only calculate the contribution for the four point function that comes from a part 
of the interaction Hamiltonian determined by the 4th order Lagrangian Hi = —L^, where L4 is the total lagrangian 
obtained from the action (|4~4")) . We should point out that the other terms that we do not consider here in the fourth 
order interaction hamiltonian are indeed important to obtain the full leading order result (see equation {75}) as was 
recently shown by fl3j . 

Of course in the end we are interested in the four point correlation function of £ and not of . At leading order in 
slow roll these two correlation functions are equal but they will differ at next-to-leading order. 



B. The trispectrum of £ 



In this subsection we calculate the relation between the trispectrum of £, on large scales, and the trispectrum of £ n 
calculated using the formalism of the previous subsection. This relation also involves lower-order correlation functions 
of £ n present in the literature. The variables £ and £„ are related up to third order by 

C(x) = C„(x) + F 2 (C„(x)) + F 3 (C„(x)), (67) 

where ^(Cn) = O(C^), ^(Cn) = C(Cn) are the terms coming from the second and third order gauge transformations 
respectively. F 2 can be found in Appendix IA1 it is 

F2(c " ) =(l + ^) c »4«" + ' s ' < 68 > 

where (3 is given in Eq. (|A6[) . In the large scale limit (super- horizon scales), we can ignore j3 as it contains gradient 
terms. F 3 was calculated in [l5| and reads 

= {A + k + ^ + t) + + 2f ) ¥ + ^ + i# + '■<« + m<-™ <»» 

where f a denotes terms that contain gradients (it can be found in [l5|). /& is the part of the third order gauge 
transformations that contains 7^ . The explicit form of ff, is to the best of our knowledge still unknown. To find out 
the explicit dependence of these terms on one would have to solve the equations of motion for 7^ , equation (I50|) . 
We do not do this in this work. We believe that these terms will vanish in the large scale limit and therefore do not 
contribute to our calculation. 

A field redefinition like C = Cn + aiCn^Cn^ + flteCn Cn j where (^■' h ' c ' d ^ denote one of ( n ,(n,(n, gives after using 
Wick's theorem a relation between both trispectrum like 

(C(x 1 )C(x 2 )C(x 3 )C(x 4 ))c - {T} + {PB} + {PPP} + 0{Pi )\ (70) 

where 



{T} = (C„(xi)Cn(x 2 )Cn(x 3 )Cn(x 4 )), 



(71) 
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{PB} 



a i 



> <~ )> + (e ) (x 1 )C„(x4))(e ) (Xl)Cn(x 2 )C„(x 3 )> 



-(e ) (x 1 )Cr l (x 3 )>(Ci Q) (x 1 )C„(x 2 )C„(x4 

+3 perm. 
+ a,2 (a — * c, 6 — * d) + 3 perm. , 



{PPP} 



^[(C^xiKi^x.))^^ 

(x 2 )Cn(x 3 )) 

+ <e ) (xi)Ci a) (x 2 ))(<C^(x 1 )Cn(x3))<e ) (x 2 )Cn(x4)) + (Ci Q) (xl)Cn(x4))<Ci b) (x 2 )Cn(x 3 )) 
+ <e ) (xi)Ci b) (x 2 ))((Ci Q) (x 1 )C„(x 3 ))(d a) (x 2 )Cn(x4)) + <C ( l Q) (xi)Cn(x 4 ))<C ( l Q) (x 2 )C„(x 3 )) 

+5 perm. 
a 2 (a ^ c, 6 — > d) + 5 perm. 

2a 1 a 2 [(Ci Q Hxi)Ci c ^x 2 ))(<CW(x 1 )C„(x 3 ))<C^(x 2 )C„(x 4 )) + (Ci b) (xi)Cn(x 4 ))(C ( l d) (x 2 )C l (x 3 )) 
+ (e ) (xi)d d) (x 2 )>((Ci b) (xi)C„(x 3 ))(Ci c) (x 2 )C„(x 4 )) + (e ) (x 1 )C„(x 4 ))(C< l c) (x 2 )C„(x 3 )) 

+ <e ) (xl)Ci C) (x 2 ))((e ) (xi)Cn(x 3 ))(e ) (x 2 )C„(x 4 )) + <e ) (xi)Cn(x4))<e ) (x 2 )C„(x 3 )) 
+ <e ) (xl)Ci d) (x 2 ))((d a) (xi)Cn(x 3 ))(Ci C) (x 2 )C„(x 4 )) + <e ) (xi)C«(x4))<e ) (x 2 )Cn(x 3 )) 

+5 perm. 



(72) 



(73) 



where "perm" means the other permutations of the preceding terms and O(P^) 4 denotes terms that are suppressed 
by successive powers of the power spectrum, (a — > c, b — > d) means terms equal to the immediately preceding terms 
with a, b replaced by c, d respectively. 

If the field redefinition contains third order pieces like C = Cn + &iCn°Vn Cn j they contribute with additional terms 



<C(X 1 )C(X 2 )C(X 3 )C(X 4 )) C = (Cn(xi)C„(x 2 )C„(x 3 )C„(x 4 )) 

+ h 



(d a) (x 1 )Cn(x 2 ))((Ci b) (x 1 )C„(x 3 ))(C( C) (x 1 )C„(x 4 )) + (Ci C) (x 1 )C„(x 3 ))(Ci b) (x 1 )C„(x 4 )) 
+ (C( b) (x 1 )C„(x 2 ))((C( a )(x 1 )C„(x 3 ))(CW(x 1 )Cn(x4)) + (Ci C) (xi)Cn(x 3 ))(Ci Q) (x 1 )C„(x 4 )) 

+ (C( C) (xi)C™(x 2 ))((e ) (xi)Cn(x 3 ))(e ) (x 1 )C„(x 4 )) + (e ) (xi)Cn(x 3 ))(e ) (xl)Cn(x 4 )) 



+3 perm. 



O(Pt) 



(74) 



To the best of our knowledge the expectation values involving operators containing derivatives of £ ra have not yet 
been calculated in the literature. However, once the mode function equation (|62p is solved, one has all the ingredients 
needed to calculate these expectation values, including the interaction hamiltonian. 



V. CALCULATION OF THE LEADING ORDER TRISPECTRUM 

In this section, we will use the formalism of the previous section and the fourth order exact interaction hamiltonian 
of subsection IIII Al to calculate the leading order trispectrum, under the assumption that the "slow-roll" parameters 
([7]|9|) are always small until the end of inflation. 
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A. The leading order trispectrum of 



To calculate the leading order trispectrum of Cn in slow roll, we need to evaluate Eq. (|66[) where Hj is read from 
the order e terms of the action (|44p . The interaction hamiltonian we get contains terms with . Fortunately it can 
be shown that to compute the leading order trispectrum we don't need to know the solution for 7^ and the knowledge 
of its equation of motion (Eq. ((59)) ) is sufficient. At this order we use the solution for the mode functions Eq. J63[) . 
The integrals in Eq. (|66p can then be performed analytically to give 



(n|C„(ki)C„(k 2 )Cn(k 3 )Cn(k4)|0) = (2 7 r)^(k 1 +k 2 + k 3 +k 4 ) 



1 



4 (ion + 3A) H2e 



Ai 



e 3 c3 n^f 
1 /'_. H 2 e 



2 6 



3A 



H 2 e 



H 2 



-An 



1 d 



2 8 



-A, 



(75) 



where the momentum dependent functions Ai are defined as 
Ai 



A' 5 ' 



^ = fc 2 fc 2 2 (k 3 -k 4 ) ^ + 3(fc 3 + h) 



I2k 3 k 4 



A 3 

(k 1 -k 2 )(k 3 -k 4 
K 



A K 2 

Z)i<7 kifcj , 3fcifc 2 /c 3 /c 4 



perm., 



A 2 



A 3 



i,.. 
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k 1 k 2 k 3 k 4 
A 3 



perm., 



(76) 



and "perm. 



refers to the 24 permutations of the four momenta. Note that the quantities of Eq. 

^ 4 k 



(|75|) are evaluated at 

the moment at which the total wave number A = Y]j—i ki exits the horizon, i.e., when Kc ssr = a*H*. This leading 
order result that comes from the 4th order Lagrangian agrees with the result of Huang and Shiu [l3| 1 that did their 
calculation in the uniform curvature gauge and using a simpler method that is only valid to calculate the leading 
order contribution for models with c 5 <l. In fact, we can compare our uniform curvature gauge result (I54[) with the 
result of Huang and Shiu [T^. We see that the last terms of the fourth, sixth and ninth lines of equation (|55|) are 
exactly the ones obtained by [l3j |. their equation (15), using the method of just expanding the field lagrangian as in 
[Til . [23] ■ For a model with a general field lagrangian these terms are the ones that give the leading order contribution 
for the trispectrum, in the small sound speed limit, equation (|75p . The contribution coming from the tensor part will 
be of next-to-leading order in this case. 

For standard kinetic term inflation, n = A = and c s = 1 and Eq. ([75]) vanishes, the leading order is then given 
by the next order in slow roll. 



B. The leading order trispectrum of £ 

It is well known that if the slow-roll conditions are satisfied until the end of inflation and we can ignore gradient terms 
then the gauge invariant curvature perturbation £ remains constant on super-horizon scales to all order in perturbation 
theory. In this subsection, we will see that this fact greatly simplifies the relation between the trispectrum of £ and 

Cn- 

In the large scales limit, Eq. (fBTj) simplifies to give 

( = Cn + a( 2 + ±Un + 0(<; 3 n ), (77) 



where a is defined as 



a=-H —. (78) 

2 2H(j)o 



1 The full leading order result for the four point function can be found in the revised version of |13[ which takes into account all the 
contributions for the fourth order interaction hamiltonian. 



14 



Using the fact that £ = on super-horizon scales and the equation resulting from a time derivative of Eq. ([77]) one 
can show that 



(79) 



This equation has a simple interpretation. The variable is not constant outside the horizon, only the gauge invariant 
quantity £ is. This is the reason why the term jjCnCn in the second order gauge transformation cannot be ignored 
when one is calculating the trispectrum of £. Substituting Eq. (f79")) in Eq. (pTj) and taking the large scale limit we 
get 



C = Cn + <n + Kn + 



(80) 



where • • • means cubic terms that contain at least one time derivative of £ ra and that will only give a contribution to 
the five point function. The variable b is defined as 



b = 



3# 2 </> 

4>Q 



ego 



eq 
3 



a 
H 



2Hfo 



~6 ' 



We shall now compare (|8H)) with the result given by the 5N formalism [28|, [2£ 
expanded in series in terms of the field perturbation as 



(81) 

31|. In the 5N approach ( is 



C = N'6<f> + -N"t 



b 2 + ^N'"6(f> 3 +O(60) 4 , 



where N is the number of e-folds and a prime denotes derivative with respect to 
the previous equation and observing that ( n = —-j-5<j> we expect 



<t>0 



1ST 



'ft 



(82) 

Now comparing Eq. (fBT)]) with 

(83) 



We verified that this is indeed the case. 

Using Wick's theorem one can now relate the connected part of the fo ur p oint correlation function of £ with the 
four point correlation function of ( n calculated in the previous section [H, [33| . This relation also involves lower order 
correlation functions of ( n , like the bispectrum (£„(xi)e n (x2)Cn(x3)) (the leading and next-to-leading order in slow 
roll bispectrum was previously calculated in 10]). The relation is 

(C( Xl )C(x 2 )C(x 3 )c(x 4 )) c = (C„(x 1 )C„(x 2 )C„(x 3 )C„(x 4 )) 

+ 2a (Cn(xi)C„(x 2 ))(Cn(xi)Cn(x3)Cn(x4)) + (Cn (*1 )(n (x 3 )) (Cn (xi )(„ (x 2 )Cn (x 4 )) 
+ (Cn(xi)Cn(x 4 ))(Cn(xi)Cn(x 2 )Cn(x 3 )) + 3 perm 

+ 4a 2 [(C„(x 1 )C„(x 2 ))(C„(x 1 )C„(x 3 ))(C„(x 2 )C„(x 4 )) 



+ (Cn(Xl)Cn(x 2 ))(Cn(Xl)Cn(x 4 ))(Cn(x 2 )Cn(x 3 )) + 5 perm 
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(Cn(xOCn(x 2 ))(Cn(xi)Cn(x3))(Cn(xi)Cn(x4)) + 3perm + 0(P, 



k > ' 



(84) 



where "perm" means the other permutations of the preceding terms and O(P^) 4 denotes terms that are suppressed 
by successive powers of the power spectrum. Now, one can easily see that at leading order in slow roll the trispectrum 
for ( and are equal, this is because the constants a and b are slow-roll suppressed. These terms will only contribute 
to the next-to-leading order corrections. 



C. The next-to-leading order corrections for the trispectrum 

In subsection IV Al we showed that for standard kinetic term inflation the leading order result, equation (|75p . 
vanishes and in fact in this case the leading order of the trispectrum of Cn is of order e~ 2 (the next-to-leading order 
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is the leading order). To obtain these leading order contributions it is easier to perform the calculation using the 
uniform curvature gauge action Eq. ([55]) - ([55]) . Eq. (|5"5|) vanishes exactly for standard kinetic term inflation. The 
action ()56|1 is exact in the slow-roll approximation but it is instructive to determine the slow-roll order of the different 
terms. One can see that the leading order contribution comes from terms of order O(e ), as pointed out in [Til ] . If we 
take P(X, <p) = X — V{4>) then the leading order (in slow roll) source of Eq. ([59]) will be of order O(e ). We therefore 
don't expect jij to be slow- roll suppressed and the actions ([57]) . ([58]) will contain unsuppressed terms of the same 
order as the leading order term of the action ([B"6]) . These tensor contributions were absent in the analysis of [ll[ and 
we have shown that they are of the same order as the fourth order action considered in [ll| , our Eq. (|56p . It is still 
an open question how these new contributions will change the trispectrum result of Seery et at. 

For the general lagrangian case, the leading order trispectrum was given in the previous subsection and in [l3| . 
Contrary to the method of 13], our method of obtaining the fourth order action ([B~f]) does not rely on any approxi- 
mation and therefore the action (|54[) is valid to all orders in slow roll and in the sound speed expansion and it can 
be used to study the next-to-leading order corrections. Depending on the momentum shape of these next-to-leading 
terms they might become big enough to be observed in the next generation of experiments. A similar argument 
applies for the next-to-leading order corrections for the bispectrum, as it was shown in [10]. For example, for DBI 
inflation, [13] showed that the leading order non-gaussianity parameters Tjvx scales like t^l ~ 0.1/cf (for a specific 
momentum configuration) and f^L ~ l/ c s- They argue that if c s ~ 0.1 then /jvl is still inside the value range allowed 
by observations but t^l ~ 10 3 could be detected with the Planck satellite CMBR experiment. Therefore, assuming 
that the slow-roll parameter e is of order e ~ 0.01 (at horizon crossing) these next-to-leading order corrections for 
the trispectrum could possible be observed with the Planck satellite. A more careful and systematic study of the 
momentum dependence of these new terms is required and it is left for future work. 



We have computed the fourth order action for scalar and second order tensor perturbations in the comoving gauge. 
Our result is exact in the slow-roll (SR) expansion but practically it is useful to study the SR suppression of the 
different terms. We were able to show that after many integrations by parts the unsuppressed terms contained in 
the previous action can be reduced to total derivatives terms plus corrections that are SR suppressed. The resulting 
action has the correct order in SR. It is suppressed by e as it should be, because in pure de Sitter space the curvature 
perturbation is a pure gauge mode. An important lesson from our work is that in order to obtain the correct SR order 
for the action, the second order tensor perturbations cannot be ignored as assumed in previous works [Txj ] and [Til ]. 
We found the explicit form of these tensor perturbations in the comoving gauge by using the gauge transformations 
from the uniform curvature gauge. Fortunately, for a general inflation model like |[T]). we showed that we do not 
need to solve the equations of motion for the tensor perturbations if we are interested in calculating the leading order 
trispectrum. However, to calculate the next-to- leading order corrections to that result, or to calculate the leading 
order trispectrum for standard kinetic term inflation, we do need to solve explicitly the equations of motion for the 
tensor perturbations. This will be left for future work. 

Using the comoving gauge action we have calculated the leading order in SR trispectrum of £■ We compared our 
result with the result of [131 ] . obtained using the uniform curvature gauge, and we found an agreement. 




For the uniform curvature gauge action, that is also exact in the SR expansion, we identified the terms that will 
contribute to the next-to-leading order corrections to the trispectrum. We pointed out that depending on the model 
and on the momentum configuration, some of these corrections might be observable with the Planck satellite. After 
taking particular limits, the previous action nicely reduces to previously know results (llj . (l3j with the caveat that 
the above mentioned works ignore tensor contributions. 

Finally we have obtained the relations between the trispectrum of £ and 5<p ( on large scales) using the third order 
gauge transformations and compared the result with the 5N formalism. 

To conclude, we have provided the necessary equations (fourth order action and the relation between £ and 8<f) to 
calculate the trispectrum for a fairly general model of inflation that are also valid for models where SR is temporarily 
interrupted, i.e., around a "step" in the inflaton's lagrangian pl|. In this case, it is impossible to apply the SN 
approach and it is required to evaluate the n-point functions numerically [l6[ (see [34| for a different approach). We 
leave this more practical application of our results for future work. 
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APPENDIX A: GAUGE TRANSFORMATIONS UP TO SECOND ORDER 



In this Appendix we will find the change of variables that one needs to perform to go from the uniform curvature 
gauge (|4"5)) to the comoving gauge (|2"Tj) . A similar result can be found in Q. In order to go from the gauge (j4"5)) 
where the field fluctuation is not zero to the gauge (f2"Tj) where <50 = we need a change of variables that satisfy 
<f>(t + T(t))+54>(t + T(t))=cj)(t). 

At first order in perturbation theory we only need to do a time reparametrization. Let t and t be the time coordinates 
in the gauges (f21j) and (|45|) respectively. The time reparametrization is t = t + T. At first order 

T = 

At second order the time reparametrization is 

T 

At this order we also need to perform a spatial reparametrization given by x % = x % + e l (x,t), where e % is of second 
order in the perturbations. The metric in the gauge (|2ip becomes 

"""WW* <S + '§ + ^ + + (* + *«> + M ■ <«> 

where iV - / is the first order shift vector in the gauge (I45[) . If the vector e l obeys the equation 

a ~ 2sh v + + = 2 /% + Mij , (A4) 

with being a transverse and traceless tensor and Shij being defined as the first four terms of Eq. (|A3[) . then the 
gauge transformation equations are given by 

HT 2 

C - HT+ — +P, 

Jij = 7« (*)+/%•• (A5) 

To obtain the quantities (3 and fiij it proves to be useful to decompose e 1 in e 1 = d l e + e\ with dit\ = 0. After a few 
mathematical manipulations of equation (|A4[) one can obtain 

= — (Shl-d- 2 d t d : >5h ij ), (A6) 



0o<50 2 50(50 



0o 20o 



(Al) 



(A2) 



Mij = ^ 2 \6hij - \s t] Sh h k - d- 2 did k 6h kj - d^d^Shki + U VJ d~ 2 d l d k 5h lk + ^d- 2 d t d 3 Sh k k + ^d-%djd l d k 5h lk j , 

(A7) 

where Shij can be written explicitly as 

11 o 2 ■ 

Sfkj = -J^di(ndj(n + — (diCndjlpl + djCndilpl) + JrCnlij, (A8) 

where is from the uniform curvature gauge and we have used the variable £ n introduced before, Eq. (|38p . As jij 
is of second order now, the term CnTij is of third order. We kept it in Eq. (|A3|) for the sake of comparison with the 
result of For £ we have 

C = Cn + k 2 n + "i^-C 2 + ^Un + P- (A9) 
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APPENDIX B: EXTRACTION OF TT PART OF A TENSOR 



Let Tij to be a given 3D symmetric tensor, as it is the case for the source of equation (|59|) . Then it can be 
decomposed into a trace part and a traceless part as 



|* y ■ /;,. (bi) 



The traceless part (5 degrees of freedom) can be written like 

Tij = D ijX + d iXj + d jXi + Xij ■ (B2) 

with Dij = didj — ^Sijd 2 , d l Xi — and d l Xij = = x\i where indices are raised by §ij. The equation d l Tj = 
\djT + %djd 2 x + 9 2 Xj can be solved using a similar method as the one we used to solve the second order momentum 
constraint previously. We then find 

X=ld- 4 d j F j , X j=d- 2 F j -d j 9 i d-^F i) (B3) 



T 

Xij = T i3 - —■),, - D^x - <>. \ i - djXi- (B4) 



where F t = d j Tj - ±<9 j; T. And 



In conclusion, given a tensor Tij, Eq. (|B4[) defines its transverse and traceless part. 

Let us see how this works at the action level for the particular case of the tensor perturbations described in the 
main text. In the action (|58p . the source for 7^ is of the form 



S = / dtd 3 xf j T i:j . (B5) 



with being quadratic in 8<j>. We can see that because 7^ is transverse and traceless we are allowed to replace 
Tj in the previous action with Xij defined in (|B4|) . If we calculate the equations of motion by varying the resulting 
action we get as a source Xij an d not simply Tj (see Eq. (|59j) ). ensuring that both sides of the equations of motion 
are transverse and traceless. 
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